home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / chpev.z / chpev
Text File  |  1996-03-14  |  4KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCHHHHPPPPEEEEVVVV((((3333FFFF))))                                                            CCCCHHHHPPPPEEEEVVVV((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CHPEV - compute all the eigenvalues and, optionally, eigenvectors of a
  10.      complex Hermitian matrix in packed storage
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CHPEV( JOBZ, UPLO, N, AP, W, Z, LDZ, WORK, RWORK, INFO )
  14.  
  15.          CHARACTER     JOBZ, UPLO
  16.  
  17.          INTEGER       INFO, LDZ, N
  18.  
  19.          REAL          RWORK( * ), W( * )
  20.  
  21.          COMPLEX       AP( * ), WORK( * ), Z( LDZ, * )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      CHPEV computes all the eigenvalues and, optionally, eigenvectors of a
  25.      complex Hermitian matrix in packed storage.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      JOBZ    (input) CHARACTER*1
  30.              = 'N':  Compute eigenvalues only;
  31.              = 'V':  Compute eigenvalues and eigenvectors.
  32.  
  33.      UPLO    (input) CHARACTER*1
  34.              = 'U':  Upper triangle of A is stored;
  35.              = 'L':  Lower triangle of A is stored.
  36.  
  37.      N       (input) INTEGER
  38.              The order of the matrix A.  N >= 0.
  39.  
  40.      AP      (input/output) COMPLEX array, dimension (N*(N+1)/2)
  41.              On entry, the upper or lower triangle of the Hermitian matrix A,
  42.              packed columnwise in a linear array.  The j-th column of A is
  43.              stored in the array AP as follows:  if UPLO = 'U', AP(i + (j-
  44.              1)*j/2) = A(i,j) for 1<=i<=j; if UPLO = 'L', AP(i + (j-1)*(2*n-
  45.              j)/2) = A(i,j) for j<=i<=n.
  46.  
  47.              On exit, AP is overwritten by values generated during the
  48.              reduction to tridiagonal form.  If UPLO = 'U', the diagonal and
  49.              first superdiagonal of the tridiagonal matrix T overwrite the
  50.              corresponding elements of A, and if UPLO = 'L', the diagonal and
  51.              first subdiagonal of T overwrite the corresponding elements of A.
  52.  
  53.      W       (output) REAL array, dimension (N)
  54.              If INFO = 0, the eigenvalues in ascending order.
  55.  
  56.      Z       (output) COMPLEX array, dimension (LDZ, N)
  57.              If JOBZ = 'V', then if INFO = 0, Z contains the orthonormal
  58.              eigenvectors of the matrix A, with the i-th column of Z holding
  59.              the eigenvector associated with W(i).  If JOBZ = 'N', then Z is
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCHHHHPPPPEEEEVVVV((((3333FFFF))))                                                            CCCCHHHHPPPPEEEEVVVV((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              not referenced.
  75.  
  76.      LDZ     (input) INTEGER
  77.              The leading dimension of the array Z.  LDZ >= 1, and if JOBZ =
  78.              'V', LDZ >= max(1,N).
  79.  
  80.      WORK    (workspace) COMPLEX array, dimension (max(1, 2*N-1))
  81.  
  82.      RWORK   (workspace) REAL array, dimension (max(1, 3*N-2))
  83.  
  84.      INFO    (output) INTEGER
  85.              = 0:  successful exit.
  86.              < 0:  if INFO = -i, the i-th argument had an illegal value.
  87.              > 0:  if INFO = i, the algorithm failed to converge; i off-
  88.              diagonal elements of an intermediate tridiagonal form did not
  89.              converge to zero.
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.